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Abstract 

Recent experiments on single crystals of YM2B2C have revealed the presence of point nodes in 
the superconducting energy gap A(k) at k = (1,0,0), (0,1,0), (-1,0,0), and (0,-1,0). In this paper 
we investigate the effects of impurity scattering on the quasiparticle spectrum in the vortex state 
of s+g-wave superconductors, which is found to be strongly modified in the presence of disorder. 
In particular, a gap in the quasiparticle energy spectrum is found to open even for infinitesimal 
impurity scattering, giving rise to exponentially activated thermodynamic response functions, such 
as the specific heat, the spin susceptibility, the superfluid density, and the nuclear spin lattice 
relaxation. Predictions derived from this study can be verified by measurements of the angular 
dependent magnetospecific heat and the magnetothermal conductivity. 

PACS numbers: 74.25. Op, 74.25.Fy, 74.70.Dd 
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1. Introduction 

The recently discovered superconductivity in the borocarbides YN12B2C and LuN^E^C 
has attracted much attention due to the rather unusual apparent structure of their gap 
function A(k). On the one hand, the presence of a substantial s-wave component has 

been observed by thermodynamic measurements of samples in which Ni was substituted by 
small amounts of Pt, causing the opening of an energy gap. 3], 0] On the other hand, there 
is also mounting experimental evidence for the presence of nodes in the superconducting 
order parameter. A consistent way to account for these seemingly contradictory 

observations is to consider the hybrid s+g-wave order parameter 0,0], given by 



A(k) = A (l - sin 4 (#) cos(40)) /2, 



(1) 



where the angles 6 and <fi indicate the direction of the quasiparticle wave vector k. A 
graphical represenation of A(k) is given in Fig. 1(a). This order parameter has rece ntly 
been shown to account for the Raman spectra measured in YM2B2C and LuNioBoC. jinl. \v\ 
Furthermore, it was observed that the corresponding quasiparticle spectrum is very sensitive 
to impurity scattering. In particular, the thermal conductivity of Y(Ni 1 _ x Pt x )2B2C for 
x=0.05 in the vortex state at T=0.8K does not exhibit any angular dependence, contrary 



to the clear cusp structure observed in the pure system. 
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FIG. 1: (a) s+g-wave order parameter with four point nodes, (b) quasiparticle density of states for 
various impurity scattering rates. 
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Theoretical studies of nonmagnetic impurity scattering in anisotropic superconductors 
indicate that this type of disorder has much more profound effects on the hybrid s+g-wave 
superconductors with point nodes than on other unconventional superconductors with line 



nodes. 



13, 



14j In particular, it was observed that the unitary limit and the weak-scattering 



Born limit give practically the same results over a wide range of parameters. Also, contrary 
to superconductors with line nodes, no low-frequency resonance is found for the strong- 
scattering case. Even more surprisingly, a disorder-induced quasiparticle energy gap uj g 
opens up even for infinitesimal scattering rates V. This dependence is well approximated by 



UJ n 



ry (i + 2T/A) 



.0 



141 ] Typical quasiparticle densities of states for the pure s+g-wave 



system and in the presence of impurity scattering are shown in Fig. 1(b). 

The opening of a disorder-induced quasiparticle energy gap offers a natural explanation 
for the reported exponentially activated temperature dependence of the thermal conductivity 
and the specific heat, 



— , — oc y— J exp(-Wg/T). 



(2) 



For the same reason, no universal heat conduction is expected for these systems as opposed 
to superconductors with line nodes. Q Because these effects of impurity scattering are 
rather singular compared with other nodal superconductors, it is necessary to reexamine 
how transport properties are affected by disorder. 0, ^] One finds that in the superclean 



limit, i.e. T <C \Jv a v c eH where v a and v c are the Fermi velocities along the crystallographic 

a- and c-directions, the thermal conductivity k zz has a similar angular dependence as in 
rn 

the pure system. |8j However, the opening of the energy gap due to the impurity scattering 
drastically reduces nodal excitations. For T /y/v a v c eH > tt/\^2 practically all of the nodal 
excitations are eliminated, and therefore the angular dependence of the thermal conductivity 
is suppressed, as it was recently observed by Kamata et al. jl^ 

The impurity scattering rate T can be accurately estimated in terms of its effect on the 
superconducting transition temperature, Q 



-In 



if) 



T c0 J l + (/ 2 ) 



1 r 

2 + 2^T 



(3) 



where ip{z) is the di-gamma function, / = sin 4 (#) cos(40) denotes the angular dependence 
of the s+g-wave gap function, (f 2 ) ~ 0.203 results from angular averaging, and T c0 is the 
critical superconducting transition temperature of the pure system. For the 5% Pt doped 



YM2B2C, the reduced critical temperature T c = 13. 5K therefore indicates a scattering rate 
T = 28. 3K, where T c q= 15.5 K was used. This implies that for the doping level x ~ 0.05 the 
corresponding disorder induced energy gap is u g ~ 10K. Since this is already a rather strong 
effect, it may be necessary to measure single crystals with x < 0.01, in order to observe the 
associated disappearance of the angular dependence in the thermal conductivity. 
2. Quasiparticle Density of States in the Vortex State 

The quasiparticle density of states and the quasiparticle scattering rate are obtained from 



the integral 



(x + l C ) ( {(x + tC ) 2 - ((T/A) + (1 - f)/2f) - 1 ' 2 ) 



(4) 



where x = (|v • q I ) / A , and v • q is the Doppler shift due to the supercurrents circulating 
around the vortices, jl 71 . 



181 ] The terms (T/A) and Cq in Eq. 4 arise from the renormalization 



of the order parameter and of the frequency due to the presence of impurity scattering, 



A 



uj=0 



iC A. 



(5) 
(6) 



For the gap symmetry of A(k) given in Eq. 1, and in a magnetic field H directed along 
, <fi) with the polar axis along the crystal c-direction, one obtains 

y/v a v c eH 



X 



2A 



1 — sin 9 cos 2 6 + \ 1 — sin 9 sin 



(7) 



For small impurity scattering T/A 1, the angular averages Eq. 4 can be performed 
analytically, giving 



/ ~ l - (x + iC Q ) I cosh" 1 + cos" 1 y - ^ I I 



y 



(8) 



where y = T/xA. Within this limit, one obtains closed expressions for the quasiparticle 
density of states g(H,T) and for the quasiparticle damping constant C A: 



g(H,T) = ^jp- = ^cos- 1 y9(l-y) + O(C ) : 



CnA 



Tx 
~2~ 



ln[|)+0(Co) 



(9) 
(10) 
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Consequently, the low-temperature specific heat, the spin susceptibility, the superfluid den- 
sity in the a-b plane, and the nuclear spin lattice relaxation rate are given by 



Cs Xs . . t \ \ 

= — = p(H,r), 



InT xn 

PSab(H.) - 3 

T^W = 1 "^ (H ' r) ' 



T -l 

1 IN 



^ 2 (H,r). 



(ii) 

(12) 
(13) 



Note that in contrast to p Sab (H) the superfluid density along the c-direction ps c (H) is not 
linear in <?(H, T). 

0.5 




FIG. 2: Dependence of the quasiparticle density of states on the azimuthal angle 4> for various 
impurity scattering strengths. The polar angle is fixed to 9 = ir/2. 

In Fig. 2, the dependence of the quasiparticle density of states <?(H, T) on the azimuthal 
angle is shown for a fixed polar angle 9 = 7r/2, and for y = 0.0, 0.4, 1.0, and 1.2, 



where yo = (2r) / (n^/v a v c eH) is a measure of the impurity scattering strength. From Eq. 
9 it follows that g(H,T) = for yo > \/2. Furthermore, it is observed in Fig. 2 that 
the ^-independent part of g(H, T) gets cut off for increasing T. For y > 1 the constant 
part is completely eliminated, and only small islands survive for 1 < y < \/2. This is 
a signature of a k-independent energy gap u g that opens up due to impurity scattering. 
Thermodynamic response functions, such as the specific heat, the spin susceptibility, the 



superfluid density, and the nuclear spin lattice relaxation rate provide a direct way to study 
this disorder-induced phenomenon. 

3. Angular Dependent Thermal Conductivity 

Following the procedure given in Ref. [3, 

the diagonal components of the thermal 

conductivity tensor are given by 



n n 21n(2/x) 



K n 2\n(2/x) 



((1 - y 2 f 2 - ^(cos- 1 y - yy/T^^j 9(1 - y), (14) 
f^j 2 (cos- 1 y' - y'y/l^f) 9(1 - y'), (15) 



where x' = a/ v a v c eH yl — sin 2 9 cos 2 0/2 A and y' = y /\ sin. 
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FIG. 3: Angular dependent thermal conductivity of a s+g-wave superconductor in a magnetic field, 
(a) zz-component, and (b) xx-component of the thermal conductivity tensor for various impurity 
scattering strengths. The polar angle is fixed to 8 = tt/2. 



The angular dependence of k zz and k xx is shown in Figs. 3(a) and (b). In the limit 
yo —>■ Eq. 14 has a similar angular dependence as given in Ref. [9|, with the exception 
of a weak hire correction in the denominator. Moreover the dependence on the impurity 
scattering strength T is similar to the one for the quasiparticle density of states described 
in the previous section. On the other hand, the structure of the angular dependence of k xx 



described in Eq. 15 is dominated by an inverted cusp centered at = it/2. Furthermore, 
we note that for small applied magnetic fields k zz increases like \fH, whereas k xx remains 
independent of H. The reason for this is that the current operator Q z vanishes at the nodal 
points in the k x — k y plane, whereas Q x is non-zero in the nodal plane. This peculiar effect 
should be explored experimentally. 

We observe that there are a number of similarities of the order parameter of these hybrid 
s+g-wave compounds with the newly discovered heavy fermion superconductor PrOs4Sbi2- 
In particular, the singlet s+g-wave gap function can describe consistently the angular de- 
pendent thermal conductivity data for this compound. 0, However, the presence of 
remnant magnetization observed in muon spin rotation measurements is suggestive of spin 



triplet pairing. 



22| More recently, Chia et al. have reported a surprisingly isotropic superfluid 



density in the B-phase of PrOs4Sbi2. 



23| It appears that none of the currently 



24, 



proposed mod- 



25 



26| Clearly, 



els for superconductivity in this material match this experimental data, 
more experimental measurements on this compound are necessary to settle this issue. 
4- Concluding Remarks 

We have seen that impurity scattering has surprisingly profound effects on s+g-wave 
superconductors such as YNi 2 B 2 C and LuNi 2 B 2 C In particular measurements of the mag- 
netothermal conductivity provide a unique window to observe the disorder induced opening 
of a gap in the quasiparticle energy through their angular dependence. Furthermore, there is 
evidence that the skuttertide PrOs 4 Sbi2 may also have a hybrid s+g-wave superconducting 
order parameter. The discovery of these point-node compounds opens up new opportunities 
to explore unconventional superconductivity. 

We thank P. Thalmeier, K. Kamata, K. Izawa, and Y. Matsuda for useful discussions on 
borocarbide superconductors, and acknowledge financial support by the National Science 
Foundation, Grant No. DMR-0089882. 
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